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We study the head-on collision of two equal-mass momentarily stationary black holes, using 
black hole perturbation theory up to second order. Compared to first-order results, this signifi- 
cantly improves agreement with numerically computed waveforms and energy. Much more impor- 
tant, second-order results correctly indicate the range of validity of perturbation theory. This use of 
second-order, to provide "error bars," makes perturbation theory a viable tool for providing bench- 
marks for numerical relativity in more generic collisions and, in some range of collision parameters, 
\Q • for supplying waveform templates for gravitational wave detection. 
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The head-on collision of two momentarily stationary black holes has become a classic problem in general relativity. 
O ., This is in part due to the fact that it is simple, when compared to more generic collisions, yet can be used as a 
D ' testbed for ideas that will later be applied in the astrophysically more relevant case of an inspiralling coalescence. 
The problem is far from academic. The near advent of gravitational wave detectors, like the LIGO or VIRGO projects, 
■ that should be able to measure gravitational radiation from colliding black holes, is a strong motivation for obtaining 
\ accurate waveform templates as soon as possible. Because one expects that most of the signal will be hidden by noise 
t-H . in a gravitational wave detector, knowing the templates accurately can mean the difference between detecting and 
not detecting gravitational waves. 

An important approach to the problem has been with numerical relativity, pioneered by Smarr § in the 70 's. 
The state of the art in these techniques yields a fairly accurate computation of radiated energies and waveforms for 
head-on (axisymmetric) coalescences || , but further development of numerical relativity will be needed to generalize 
computations to cases with less symmetry. Recently [^,^| , perturbative techniques applied to this problem were found 
\Q I to give remarkably good results when the holes were initially close. By treating the collision of close black holes as 
ON ■ the dynamics of a single distorted black hole, one can use the well understood machinery of black hole perturbation 
| theory to treat the problem. Better yet, generalization of these techniques to collisions without axisymmetry entails 
Qi only minor modifications. 

A key missing element in perturbative close-approximation calculations is that one cannot tell a priori over what 
(3JT), range it is valid. (How close is "close enough"?) The main point of this article is to show that this missing element 
is provided by second-order perturbation theory, which has recently been developed ||, largely for this purpose. A 
higher-order calculation is characterized by a number of interesting technical issues (such as gauge choices of different 
order) , and by considerable complexity. Here we want only briefly to describe the nature of the computation, and to 
argue the importance of the first results. We initially specialize to the details (e.g., even-parity only) of a head-on 
collision of two non-spinning equal-mass holes starting with Misner |7| initial data. This is done for definiteness in 
the description, and because this is the case which can most thoroughly be compared with the numerical results. 

The perturbation calculation starts with an expansion of the metric, in terms of a perturbation parameter e, in the 
form g^y — <?^ hw + eh$ + e 2 h^ + 0(e 2 ), where <?^ hw is the usual Schwarzschild metric, for a spacetime of mass M. 
Work on perturbation theory, more than twenty years ago, showed that for even-parity perturbations all the physical 
information in each ^-pole of the first-order perturbations h^J can be encoded into a single "Zerilli function" ip^ that 
is a combination of the h$ and their derivatives |^||. This function satisfies the Zerilli equation, a simple, linear 
wave-like equation. If we construct essentially the same combination of h^J to form a second-order Zerilli function 
w , the Einstein vacuum equations, to second order in e, guarantee that ip( 2 ^ will satisfy a wave equation 

- d 2 4 2) + d^4 2) + Veir)^ = S e , (1) 

which differs from the usual Zerilli equation only in the presence of a nonzero "source" term on the right fllCfl . In this 
equation, r* = r + 2M\n(r/2M — 1) and Ve(r) is given, e.g., in H||. The new feature for second-order, the source 

term Se, is a quadratic function of ip£ . In the case of the head-on collision, though not more generally, it turns out 
that only the 1 = 2 first-order terms contribute to the quadrupole source S2- For this case the source is explicitly 
given in M. 
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There are two ways in which il^p is not exactly the same combination of the hffi, as ip^ is of the h^u- First, for 
technical reasons ||, is actually equivalent to the time derivative of the second-order Zerilli function. Second, 
V^ 2 ** is not unique; one can redefine it up to terms quadratic in the first-order perturbations. The choice made in |J 

(2) C2) 

ensures that the source term is such that it yields a well behaved ip 2 m the radiation zone. (Here ^ corresponds 
to x hi the notation of ||.) 

For the head-on collision problem, the perturbation expansion parameter e is a measure of the initial separation of 
the holes. The particular choice we make is the parameter denoted k 2 in Ref. j| that occurs rather naturally in the 

description of the problem. It is straightforward to expand the Misner initial data in this parameter and to find h^} 

(2) ...... . 

and Kpv 1 from which we derive the initial conditions for the dominant 1 = 2 perturbations, 
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_ ,( 2 ), 16384 M 6 [5(r - M)^r - 2M + 19^r(r - 2M)] 
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(The surprising nonzero time derivative in (0) arises partly because ip^ is actually the time derivative of the Second- 
ly 

order Zerilli function, and partly due to the quadratic terms that were added to ip^ to make it well behaved in 

the radiation zone [pT|.) The solution for , an d therefore the source term for ([!]), is known numerically from the 
solution to the first-order problem jij. This source and the initial data of (||) and (|J) complete the specification of 
the problem defined by (|l|). 

With the solution of that problem for 0>rj , horn a simple finite differencing scheme (and with tf>^ from the 
first-order solution), the radiated power is found from 
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and the energy radiated is the time integral of this expression. The expression in brackets {} is proportional to our 
second-order radiation wavefunction, and omits yet-higher-order corrections. In applying (^), it would be inconsistent 
to keep the terms of order e 4 , which are of higher order than omitted terms, so we compute energy only to order e 3 . 

Figure 1 shows the energy results (to order e 3 ) as a function of the Misner parameter /i . This figure compares 
these second-order results with first-order results and with those of numerical relativity, and illustrates the point we 
wish to emphasize. If we had only the first-order calculations, we would know that the predictions at very large 
separation parameter were unreasonable, but we would not know how small no must be for the results to be 
reasonably accurate. By comparing the first-order prediction with that of the second-order computation, we can 
infer that perturbation theory is applicable up to /iq of order 1.8 or so. (The actual cutoff would depend on what 
level of accuracy one demands.) In the case of the head-on collision from Misner data, this prediction is subject 
to verification by comparison with the results of numerical relativity. Figure 1 clearly shows that the prediction is 
correct; perturbation theory is valid (i.e., agrees with numerical relativity results) up to around (1q = 1.8. Second- 
order calculations do in this case provide "error bars" for perturbation calculations. There is nothing specific about 
the head-on problem that favors this outcome, and it is reasonable to infer that that second-order results will provide 
the same "error bars" in generic coalescences. 

A mild caveat hangs over the application of this method to other cases. There is no unique division of perturbations 
into orders ]l2| ]. If, for example, one knew the exact radiated energy for the head-on collision, for all values of /j,q, 
one could define a new expansion parameter e that depends on /io in such a way that first-order perturbation theory 
would be exact! More worrisome, in any perturbation problem one could define an e for which the first-order energy 
is not exact, but for which the second-order perturbations are zero. The "error bars" provided by these second-order 
perturbations would be wildly misleading Jl3[ |. Such cases keep our method of error-bar determination from being 
rigorous, but are a problem more in principle than in practice. Anomalous behavior of an order of perturbations arises 
only in atypical circumstances, usually when one contrives to arrive at just such an atypical result. For a "generic" 
perturbation parameter one expects the second-order error-bar method to be just as useful as in Fig. 1. 
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A secondary advantage of second-order computation is that it can give improved accuracy of results. The difference 
between first-order results and second-order-correct results is, of course, significant only for cases in which perturbation 
theory is marginally applicable. But it is just these cases for which the radiation is strongest and which will typically 
be of greatest interest for comparison with numerical relativity, or for considerations of detector design. In Fig. 2 
we show a comparison of waveforms for the marginal case of /iq = 2, corresponding to an initial proper distance 
L = 3.3Af between the throats. At this separation there is no single all-encompassing horizon surrounding both 
throats, so the assumptions of the close limit approximation are violated. As expected, the first-order result shows 
significant disagreement with the numerical relativity result, but the result of our second-order computation is a 
waveform in remarkably good agreement with the numerical relativity result. 

The success is not only remarkable, it is at first somewhat puzzling. The waveform for second-order theory has an 
amplitude that is at most around 10% larger than the numerical relativity waveform, implying that the second-order 
result for radiated energy should be around 20% larger than that from numerical relativity. In Fig. 1, however, the 
second-order energy is about half the result of numerical relativity at /io = 2. The explanation of this brings out an 
interesting point about the computations. If we take the second-order waveform of Fig. 2, square it, multiply it by 
the correct factor and integrate, we calculate an energy which is, as it must be, about 20% larger than that from 
numerical relativity (i.e., the result of the same operations on the numerical relativity waveform of Fig. 2). But this 
procedure for computing second-order energy includes the formally disallowed e 4 terms in (^). The success of this 
formally inconsistent procedure suggests that the omitted terms in (^|) are small, so that the effect of omitting them 
is unimportant. We must keep in mind that for a marginal case like fi = 2 the expansion parameter is of order unity, 
so that "order in e" is not necessarily the critical factor in the importance of a term. The omitted terms in (^) turn 
out to be the result of a nonlinear mixing of £ — 2 and £ — 4 contributions. The £ = 4 contribution for the problem 
is small, quite aside from the fact that it is higher order in e. [See Fig. 12 of Ref. ||; at /i = 2, the £ — 4 energy is 
three orders of magnitude smaller than the £ = 2 energy] Our higher-order waveforms are more accurate than are 
implied by a comparison of first-order and second-order-correct energies, but this fact may be peculiar to the specific 
case of initial Misner data, and not applicable to more general black hole collisions. 

It is the study of more general collisions to which this work is directed. We have shown that we have a practical 
and reliable (though not rigorous) index of the accuracy of perturbation theory, and hence the basis for a consistent 
investigation of collisions that are of interest as sources of detectable gravitational waves. With estimates of accuracy 
established by higher-order calculations, perturbation theory will give waveforms and energies that will check numerical 
relativity codes; that can help in exploring the parameter space of collisions for those cases of interest for greater 
numerical work; that can in some cases provide guidelines for development of details of detection schemes. Work is 
already underway on the application of this technique to the head-on collision of initially moving and spinning holes, 
and has started on the modifications that are needed when the final hole is rapidly rotating. 
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FIG. 1. The total radiated energy in a black hole collision as a fraction of the spacetime mass M for different values of the 
initial separation. The upper scale is in units of L/m, where m is the mass of each individual hole (inferred from the minimal 
area of the throat for each hole); the bottom scale is in the parameter fio that appears in the Misner data. The dotted curve is 
the first-order result and the solid curve the result up to second order. For hq ps 2.2 the second-order corrections are as large 
as the first order result. The numerical results (with error bars) correspond to the NCSA simulations of Ref. [3]. 
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FIG. 2. Waveforms for /io = 2. Shown are the waveform for time derivative of the first-order Zerilli function and that 
waveform with the second-order correction added. This is compared with the equivalent waveform from numerical relativity. 
We see a significant improvement due to the use of second-order perturbation theory, even for this case where the separation 
of the two black holes is L/M » 3.3. The inset shows late oscillations. The period of the numerical relativity result for these 
oscillations becomes larger than the quasinormal period and is thought to be a numerical artifact [5] . The second-order correct 
curve then is probably more accurate in this region. The full numerical data sets have higher resolution than the one depicted, 
we have omitted some data points for visual clarity. 
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